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Hydrogenated nanographite can display spontaneous magnetism. Recently we proposed that 
hydrogenation of nanographite is able to induce finite magnetization. We have performed theoretical 
investigation of a graphene ribbon in which each carbon is bonded to two hydrogen atoms at one edge 
and to a single hydrogen atom at another edge. Application of the local-spin-density approximation 
to the calculation of the electronic band-structure of the ribbon shows appearance of a spin-polarized 
flat band at the Fermi energy. Producing different numbers of mono-hydrogenated carbons and di- 
hydrogenated carbons can create magnetic moments in nanographite. 

PACS numbers: 75.75.+a, 75.50.Xx, 73.20.At, 73.22.-f 
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Recent reports on magnetic carbon have stimulated 
renewed interest in carbon systems as possible new mag- 
netic materials made exclusively of light elements. These 
investigations include experiments suggesting ferromag- 
netic behavior in polymerized fullerenes made using 
photo-assisted oxidation H or high-pressure treatment. 
A nanometer-scale graphite called nanographite may 
yield another magnetic carbon system. Some ex- 
periments have demonstrated the existence of anoma- 
lous magnetic behavior unexpected in diamagnetic bulk 
graphite. §|| 

Several theoretical investigations of nanographite 
have been carried out prior to or parallel to these 
experiments. §, @, g, |, [□] |f, || |]j An important 
finding was that 7r electrons in nanometer-scale graphite 
are strongly affected by structure of edges. A typical 
edge is the zigzag edge. The zigzag edge was found in 
STM images of nanometer scale graphene, [[l5j |l6| al- 
though the armchair edge has lower energy than the 
zigzag edge. |l7], [ll| Fujita and coworkers suggested that 
the 7T electrons on a mono-hydrogenated zigzag edge may 
create a ferrimagnetic spin structure on the edge.J^] Mag- 
netism with localized spin moments is possible due to the 
existence of non-bonding localized states at the zigzag 
edge. These states are called edge states and are highly 
degenerate at the Fermi level and may be spin polar- 
ized. Such an edge state does not appear on an armchair 
edge. Klein has studied another edge of graphene strips, 
considering a zigzag edge with an atomic site having a 
7r-orbital on the edge.[fl9| po| Examining a tight-binding 
model of n electrons, he found degenerate surface states 
on the edge, which is referred to here as a bearded edge. 

The total spin moment of a graphene ribbon, however, 
should be zero, if both of the two edges of the ribbon are 
zigzag edges, or if both of them are bearded edges. This is 
because, the structure of the ir network becomes a bipar- 
tite lattice having the same number of sublattice sites. In 
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such cases, an argument using the Hubbard model results 
in zero spin moment, following the Lieb theorem. || |2l[] 
Naively, we may say that local magnetic moments at the 
two edges are coupled anti-ferromagnetically. Q 

In this paper, we show that a hydrogenated graphene 
ribbon can have a finite total magnetic moment. A key 
point is the structure of the two edges of the ribbon. 
In our ribbon, one of the edges is composed of mono- 
hydrogenated carbon atoms and another is made of di- 
hydrogenated ones. By a first-principles electronic-state 
calculation, we will show that the system is stable and 
has a fully spin-polarized flat band. A finite spontaneous 
magnetization appears in the graphene ribbon. 

We present here our idea to find magnetic 
nanographite. If a graphene ribbon has a bearded 
edge on one side and a zigzag edge on the other side, 
the 7r network becomes a bipartite lattice with different 
numbers of sublattice sites. [^2], |23| Then, a flat band 
appears at the Fermi level. The band is fully spin 
polarized if inter-electron interactions e.g. the Hubbard 
repulsion U is considered. pi ] In addition, the flat band 
consists of edge states, so that the magnetic moments 
are localized at both edges and the distribution of 
moments shows exponential decay in the tt network. p3| 
However, to design a realistic material, we also have to 
specify the atomic configuration of the bearded edge. 

We found that one solution is a zigzag edge whose 
end carbon atoms are di-hydrogenated.(24| In order to 
present an example of magnetic nanographite, we will 
show the electronic band structure of the ribbon which 
has a completely spin-polarized flat band at the Fermi en- 
ergy. The magnetism will be explained using a Hubbard 
model. 

To investigate the stability and magnetic structure of 
the nanographite, we used an electronic-state calculation 
with a local spin-density approximation (LSDA) based 
on the density functional theory. p5j, p(| An LSDA func- 
tional given by Perdew and Wang |27|| was utilized as the 
exchange-correlation energy functional. The ultra-soft 
pseudopotential generated using the Vanderbilt strategy 
was adopted, p^] The valence wave functions were ex- 
panded in a plane- wave (PW) basis set. To achieve con- 
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vergence in the total energy, we used an energy cut-off E c 
of 49Ry for structural optimization in a fixed unit cell. 
To reduce computation time, we adopted E c — 25Ry to 
optimize the size of the unit cell before the final opti- 
mization of inner coordinates. Structural optimization 
was performed until each component of the interatomic 
force became less than 8 x 1CP 5 H/a.u. The calculations 
were done using a computer program called the Tokyo 
ab-initio program package (TAPP).[p9[ 

To check the accuracy of our calculation, we obtained 
the structural parameters of a hydrogen molecule, a CH4 
molecule, a C2H4 molecule, diamond and graphite. The 
bond length of each molecule and the bond angles of CH4 
and C2H4 were obtained with relative errors less than 
3.3%. The vibrational energy hui e of H2 was 1.92 x 1CP 2 
[H], which is about 4% smaller than 2.00 x 10" 2 [H] ob- 
served. For cubic diamond, the lattice constant, the bulk 
modulus and its pressure derivative were 6.677 a.u., 4.66 
Mbar and 3.68. For hexagonal graphite, the lattice con- 
stants were a = 4.616 a.u. and c = 12.67 a.u. These 
values are in reasonable agreement with the values de- 
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Band structures of the solids also reproduced previous 
results. |§ 

Consider the graphene ribbon shown in Fig. [l]. In this 
structure, each dark circle represents the position of a 
carbon atom. Solid lines represent a bonds. Each open 
circle is a hydrogen atom. For convenience, we number 
atoms and name them CW • • • and H« ■ • • H< 3 ' 

in a unit cell. Here, n is a positive integer. On the left 
edge, each carbon atom is bonded to a hydrogen atom 
and has a tt orbital. On the right edge, each carbon atom 
has bonds to two hydrogen atoms and two neighboring 
carbon atoms. If we consider that the di-hidrogenated 
carbon atom C^ 2 "* 1 forms sp 3 bonds instead of sp 2 bonds, 
C (2n) loses a 7T electron. Since there remain In— 1 carbon 
atoms having a ir electron per each in a unit cell, the 7r 
network might become a bipartite lattice with different 
numbers of sublattice sites. 

We can consider ribbons with various widths. How- 
ever, we know that the optimum width is given by 
n ~ 5 ~ 10 for nanographite ribbons with zigzag edges. |7j 
In any case, the essential properties should be the same 
irrespective of the width. Thus we first consider a ribbon 
with n = 5 and discuss dependence of magnetism on the 
width of the graphene ribbon at the end of the paper. 

To examine our idea, we have performed an LSDA cal- 
culation. Below, we assume a periodic boundary condi- 
tion. Since the structure is one-dimensional, we are al- 
lowed to use k points along the y axis as (0, k y , 0) for the 
structural optimization. We used 6 k points uniformly 
distributed in the Brillouin zone for most simulations. 
We determined the equilibrium lattice constant of the 
ribbon, which became 4.674 a.u. For the calculation to 
optimize the inner parameters, the size of the cell was 
taken to be 40.27 x 4.647 x 18.60 (a.u.) 3 . The spacing 
between neighboring ribbons became 19.5 a.u. in the x 
direction and 15.3 a.u. in the z direction. The structure 




FIG. 1: A structural model of a magnetic graphene ribbon. 
Dark circles represent carbon atoms and open circles denote 
hydrogen atoms. Dashed lines represent a section of a unit 
cell in the xy-pl&ne. The unit cell contains 2n carbon atoms 
and 3 hydrogen atoms. Here a case with n — 5 is shown. Two 
sets of 3-dimensional coordinate axes, (xyz) and (x'y'z'), are 
displayed. The a;-axis and the t/-axis are in a graphene plane. 
Three axes x', y' and z' point to the (011), (0-11), (100) 
directions, respectively, in the {xyz) coordinate. 



was optimized by a conjugate gradient method. The ob- 
tained results did not change when we used 12 k points 
instead. The bond length of C-C ranged from about 2.64 
~ 2.70 a.u. for neighboring carbon pairs except for 2.59 
a.u. between and C< 9 > and 2.76 a.u. between C*- 9 - 1 
and C^ 10 ). Bond angles H( 2 )-C( 10 )-H< 3 ) and C^-C^ 10 )- 
are 100.7° and 114.8°. These structure parameters 
suggest that the end carbon C^ 10 ^ forms sp 3 bonds. 

An electronic band structure of the graphene ribbon 
is shown in Fig. |[ In the first Brillouin zone, the T-Y 
line corresponds to a one-dimensional axis on which a 
k vector is parallel to the ribbon. A flat band appears 
at the Fermi level. The band is almost dispersionless in 
the whole Brillouin zone and is completely spin polar- 
ized. That is, only the spin-up branch is filled and the 
spin-down branch is empty. Splitting between polarized 
flat bands for up-electrons and down-electrons is approx- 
imately 0.5 ~ 0.6 eV and the gap between them is ~ 0.2 
eV. These values give a rough energy scale of stiffness for 
the magnetic ground state. The total spin Stot becomes 
1/2 times number of unit cells N c . Stability of the polar- 
ized flat band was confirmed by using a unit cell which 
was doubled in the y direction. 

The flat band is the center band of 9 well-characterized 
7r-bands. See band dispersion around the Y point in Fig. 

The dispersion relation of these 7r bands is roughly 
the same as those obtained for a simple tight binding 
model with 9 7r-orbitals on a graphene ribbon. In Fig. |^, 
all branches for up-electrons and down-electrons are split 
by spin polarization. This result implies that a magnetic 
moment appears on every atom. We show obtained spin 
density in Fig. H. A ferrimagnetic spin structure is clearly 
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FIG. 2: The electronic band structure of a hydrogenated 
graphene ribbon. Solid lines and dashed lines represent the 
dispersion relation of bands for up electrons and down elec- 
trons, respectively. At the Fermi level Ef — 0, a dispersion- 
less band appears. The band is completely spin polarized. 
Above and blow the flat band, there are 8 7r-bands, which are 
from the 18th band to the 21th band and from the 23th band 
to the 26 th band on the K-Y line. The 17th band on the 
K-Y line is composed mainly of Is orbitals of two hydrogen 
atoms H (2) and H (3) in Fig. 1. 




FIG. 3: The optimized structure of a magnetic graphene rib- 
bon. Yellow spheres and white spheres represent the positions 
of carbon atoms and hydrogen atoms, respectively. Spin den- 
sity on the ribbon is displayed by isosurfaces. Red surfaces 
and blue surfaces represent spin-up density and spin-down 
density. The isolevel is set to 1/4 of the maximum amplitude 
of the spin-up density. 



seen. A local spin moment exists on every carbon atom. 
Interestingly, finite spin moments appear on two hydro- 
gen atoms which have bondings with C^ 10 ^. The reason 
will be explained below. Thus, the graphene ribbon is 
confirmed to be magnetic in the LSDA model. We call 
this structure magnetic nanographite. 

We found a curious characteristic of the wave functions 
of the flat band. We surmise that only 9 carbon atoms in 
the unit cell possess 7r orbitals which contribute to the 7r 
bands. This is because sp 3 bonds are created at the di- 
hydrogenated carbon atom C^ 10 '. Wave functions on the 
flat band were non-bonding edge states with finite am- 
plitude on only odd-numbered carbon atoms, but not on 
even-numbered ones. An interesting point is that a finite 
amplitude was found on two hydrogen atoms H^ 2 ) and 
H(3). Therefore, we have to consider Is orbitals on these 
hydrogen atoms in order to analyze the spin-polarized 



flat band. 

Let us introduce two sets of orthogonal coordinates, 
(xyz) and (x'y'z 1 ) as defined in Fig. |l|. Four sp 3 orbitals 
on C^ 10 ) may be represented as 4>\ — 4> s + <fi x > + </y + <f> x > , 

<t>2 = 4>s - <f>x' - 4>y' + <t>x>, 03 = 4>s + 4>x> ~ 4>y> ~ <f>z' , 

and 04 = 4>a ~ 4>x' + <t>y> — 4>z' ■ Here, <p s is a 2s orbital 
and 4> m [m = x',y',z') are three 2p orbitals on C' 10 -'. 
4>i and (j)2 point to and H^ 3 ', respectively. We take 
linear combinations of <pi and <f>2 as well as Is orbitals <p\ 
and ip 2 on H^ 2 ) and H^ 3 ) and define 4>± — (cj>\ ± (/) 2 )/v / 2 
and (p± = (<pi ± </?2)/v / 2- Note that </>_ = </> 2 , i.e. a 
7T orbital on C*- 10 ^. By symmetry, c/>_ and (p>- hybridize 
with the 7T orbitals of the host 7r system of the graphene 
sheet, but not with a orbitals in the honeycomb lattice. 
Therefore, we must utilize </>_ and (p- as well as the it 
orbitals on CW (i = 1, • ■ ■ 9) to analyze the tt network in 
the magnetic nanographite. 

If we adopt the above representation, we obtain an 
effective model describing the ribbon. The model consists 
of 11 local orbitals in a unit cell. Here we ignore overlap 
integrals between orbitals for simplicity. Introducing the 
operator representation Ci tG for a 7r electron on a carbon 
atom and d; i<T for an electron in (p~- in the Z-th unit cell, 
we obtain an effective Hamiltonian. 

H = -J2J2Uj(cl a c 3 ^ + H.c.)+J2Un zA n hL 

cr I 

+ £ ^2^2 m i,T (i) 

Here, an index i specifying a carbon atom runs over all 
carbon sites. Two kinds of number operators are defined 
as n, i<T = c\ a Ci tC , and m/ i<T = d\ a di^. Bond connec- 
tions given by non-zero £y and i! i l are the same as those 
in a grap hene ribbon with a zigzag edge and a bearded 
edge.[g3| We set the origin of the energy at the orbital 
energy of the 2p orbital e^-p on carbon. E\ s denotes the 
energy of a Is electron on a hydrogen atom and e is given 
by e = E\ s — £2p- Since e < 0, a band mainly composed of 
(p- is separated in energy from other tt bands. This prop- 
erty is also seen in the LSDA band structure, where the 
17th band on the K-Y line in Fig. ||is the corresponding 
band. 

Consider a periodic boundary condition for a finite 
system with iV c unit cells. If e = and if U' = U, 
Lieb's theorem is applicable and the ground state be- 
comes a ferrimagnetic state with a total magnetization 
Stot = \Na — Nb\/2. Here the difference in sublattice 
sites |A0l — Nb\ equals N c . This magnetic ground state 
is thought to be stable against perturbation like e < 0, 
since 1) there is a finite gap between the flat band and 
other 7r bands both in the tight-binding description and 
in the LSDA result in Fig. 2) we have confirmed that a 
small difference of ei s and S2 P does not break the flatness 
of the center band in the tight-binding model. Thus the 
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magnetism found in an LSDA calculation can be under- 
stood as the appearance of ferrimagnetism in a flat-band 
Hubbard model. 

The formation energy of a di-hydrogenated edge has 
been evaluated by comparing the total energy within 
LSDA. We denote Eh 2 , Eqx and Eg2 as the energy 
of a hydrogen molecule, that of a graphene ribbon with 
mono-hydrogenated zigzag edges and that of a magnetic 
nanographite ribbon, respectively. The latter two are 
energy per a unit cell including 2n carbon atoms. Our 
calculation shows that E G2 — {Eqi +E H2 /2) ~ 0.3eV for 
n = 5. Thus, the hydrogenation of nanographite with 
mono-hydrogenated zigzag edges is confirmed to be an 
endothermic reaction. Utilization of proper catalysts in 
catalytic reduction may be important in the synthesis of 
magnetic nanographite. 

To see how our conclusion is affected by changing the 
width of the ribbon, we performed the LSDA calculation 
for ribbons with a width from n = 2 to n = 4. A ferrimag- 
netic structure was found even in these narrow ribbons. 
These ribbons do not show the Peierls distortion, which is 
similar to conclusions on the mono-hydrogenated zigzag 
ribbons. |, [u| || 

We have shown that a spin-polarized flat band appears 
in the LSDA band structure of a graphene ribbon. The 



origin of the magnetism is explained as a realization of 
the flat-band ferromagnetism. An orbital (p_ on two hy- 
drogen atoms H^ 2 ) and H^ 3 ) is an anti-bonding orbital 
and has symmetry p z . Thus, hydrogenation of a car- 
bon atom is identified as the addition of a 7r orbital at 
a zigzag edge. The structure may be regarded as a real- 
ization of the bearded edge. We can expect magnetism 
for any hydrogenated ribbon with \Na — Nb\ > 0, even 
if di- hydrogenation of the zigzag edge is not perfect. The 
same effect is also expected with fluorination. [^2) In this 
report we have considered a ribbon with a flat graphene 
structure, but the method is also applicable to nanotubcs 
with zigzag edges. Thus it appears that, a magnetic nan- 
otube may be designed. p4| 
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